Pseudopotential model of ultracold atomic collisions 
in quasi-one- and two-dimensional traps 
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We describe a model for s-wave collisions between ground state atoms in optical lattices, consid- 
ering especially the limits of quasi-one and two dimensional axisymmetric harmonic confinement. 
When the atomic interactions are modelled by an s-wave Fermi-pseudopotential, the relative motion 
energy eigenvalues can easily be obtained. The results show that except for a bound state, the trap 
eigenvalues are consistent with one- and two- dimensional scattering with renormalized scattering 
amplitudes. For absolute scattering lengths large compared with the tightest trap width, our model 
predicts a novel bound state of low energy and nearly-isotropic wavefunction extending on the order 
of the tightest trap width. 

PACS numbers: 32.80.Pj, 32.80Lg, 34.50.-s 



I. INTRODUCTION 

Unprecedented control has been gained over atomic 
collisions through cooling to nanoKelvin temperatures 
and selection of internal hypcrfinc states. Further con- 
trol of these systems is now being gained by manipulating 
the atoms with external fields, including magnetic fields 
and optical lattices. An opportune example where both 
collisions and the external trapping potential are essen- 
tial is the superfluid-Mott insulator transition in an op- 
tical lattice. This transition was recently demonstrated 
by beginning with an atomic Bose-Einstein condensate 
and adiabatically turning on an optical lattice 0. The 
final Mott insulator state has a fixed number of atoms 
per lattice site. Such a system represents an ideal en- 
semble for measuring scattering properties of the atoms. 
The free-space scattering amplitudes determine the en- 
ergy eigenvalues and loss rates of the system and vice 
versa. 

There are many current and proposed applications of 
cold atoms in optical lattices. One is high-precision mea- 
surement of atomic potentials, by determining positions 
of weakly-bound molecular states for example 0|- An- 
other is the implementation of quantum logic gates with 
neutral atoms. Several proposed implementations com- 
bine optical lattices with internal-state-dependent cold 
collisions 0, l^- One can also consider the combina- 
tion of an optical lattice along one or two directions 
and a weak dipole trap in the remaining directions. 
Such highly anisotropic trapping configurations were al- 
ready used in experiments on Bose-Einstein condensate 
number-squeezing Q and controlled loading |]^, as well 
as looking for dimensional effects on the condensate den- 
sity J]. Optical '8^ or magnetic |9j waveguides have been 
used for guiding cold atoms. These quasi-one- and two- 
dimensional configurations can be combined with a tun- 
able scattering length, such as from magnetic Feshbach 
resonance, to lead to new physical regimes. Workers on 
Bose-Einstein condensate experiments are observing in- 
creased phase-fluctuation in the crossover to one dimen- 



sion [loL[l]|. Theoretically, the (one-dimensional) Tonks 
gas regime can occur at low density when bosons are 
tightly confined along two directions and weakly con- 
fined along the third [l^ IT^ . One- half anyon statis- 
tics and the fractional quantum Hall effect may be pos- 
sible for bosons confined tightly along one direction and 
weakly along the other two Another proposal is 

the Kosterlitz-Thouless transition for interacting bosons 
in two-dimensions 15]. For all these applications, one 
first needs to understand the relation between cold colli- 
sions in free space and in a trap with arbitrary geometry. 
Then the many-body physics can be treated on the ba- 
sis of effective low-dimensional interactions. In this ar- 
ticle, we show how to compute the eigenvalues from the 
free-space scattering data and trap frequencies for arbi- 
trary axisymmetric harmonic traps, emphasizing the one- 
and two-dimensional regimes. We also describe a novel 
bound state induced by the trap in both low-dimensional 
regimes when the scattering length is large in magnitude. 

Previous theoretical work has included exact solu- 
tions of collisions for special interatomic potentials in an 
isotropic harmonic trap IT^ . and comparison with 
results including realistic ground-state interatomic po- 
tentials. We have shown how an effective-scattering- 
length combined with a Fermi pseudopotential can be 
used to calculate the energy eigenvalues for collisions in 
an isotropic harmonic trap 0, E|- We compared our 
model with full interaction potential results for both a 
single-channel collision and a multi-channel collision with 
a magnetically tunable Feshbach resonance. In both 
cases, the model can accurately treat tight traps, as long 
as the trap size is larger than the van der Waals scale 
length. On the strength of this evidence we propose to 
also apply the effective-scattering length model to treat 
ground state collisions in axisymmetric harmonic traps. 

The effect of confinement along only one or two di- 
mensions has also been considered theoretically by oth- 
ers. For one dimension, the scattering could be rep- 
resented by a one-dimensional delta function provided 
the coupling strength was renormalized by the confin- 
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ing trap [T5| . Quasi- two-dimensional scattering solutions 
were found for two dimensions, where the renormaliza- 
tion is of the two-dimensional scattering length ■ We 
will consider the applicability of these two results in the 
case of very prolate and oblate traps respectively. 

The outline of our paper is as follows. In the next 
section, we state the problem of atoms colliding in an 
axisymmetric trap and the regime of applicability of our 
effective-scattering-length model. In Sec. IIIII we show 
how this problem can be solved in a particular basis. 
Sec. IIVI contains the results of the calculation for the 
one-dimensional regime and a comparison with a one- 
dimensional scattering theory; Sec. Ivl contains the same 
but for the two-dimensional regime. In Sec. lVII we discuss 
a novel trap-induced bound state which appears in both 
the one- and two- dimensional regimes for large positive 
or negative scattering length. We give an example of 
magnetically tunable Feshbach resonance for Na atoms 
in quasi-one-dimension in Sec. lVIll We conclude the pa- 
per in Sec. IVIIII In the appendix we derive the matrix 
elements of the anisotropic potential required for the cal- 
culation. 



Virit(r), and has the usual harmonic oscillator solutions. 
The relative coordinate Hamiltonian is given in spherical 
coordinates by 



H 
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2^r^ dr 
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dr 



V{r,e) + V,nt{r) (4) 



where r = ri — r2 and r — |r| . The interatomic orbital an- 
gular momentum operatorL gives the partial wave quan- 
tum numbers / = 0, 1, 2, ... for s-, p-, d-, ... waves. The 
potential due to the trap written in terms of the spherical 
harmonic Y2Q is 



V{r,e) = ^nuj^r^ 
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where the mean-square trap frequency appears as 



■UJt 



(5) 



(6) 



and 



II. TWO ATOMS COLLIDING IN AN 
ANISOTROPIC HARMONIC TRAP 

We assume ultracold atoms are trapped in an optical 
lattice detuned far off resonance. Specific optical lattice 
potentials for different angular and polarization configu- 
rations are calculated in 20] . For our purposes we assume 
that two atoms in specific internal sublevels remain near 
a local minimum of the potential. With the assumption 
of local azimuthal symmetry about an axis through a 
potential minimum, we approximate the anisotropic po- 
tential near a particular site by 

Vt{v,) = \m {^\{x] + vf) + iol??^ , (1) 

where Vj is the positions of atom j = 1 or 2, m is the 
atomic mass, and a;_L and Wz are trapping frequencies. 
We define the trap anisotropy 



A = — 



(2) 



so that one-dimensional physics is approached for the 
oblate trap A ^ 1 (" cigar" ) and two-dimensional physics 
for the prolate trap A 3> 1(" pancake"). The length scales 
associated with the transverse and longitudinal trap di- 
rections are 



di = 




— ^ (3) 



where /.t = m/2 is the reduced mass of the atom pair. 

As in the case of an isotropic harmonic trap, the center 
of mass and relative motion are separable. The center of 
mass motion is independent of the interatomic potential. 



■2ujI 



A^ 



(7) 



The term proportional to A in Eq. ^ defines the 
anisotropic part H^^") of the Hamiltonian H = H^^^ + 
H^^\ The length scale associated with the mean-square 
frequency is defined to be 



(8) 



The interatomic Vi-nt{r) potential for two ground state 
atoms approaches —C^/r^ at large internuclear separa- 
tion r. The associated van der Waals length scale is 
xq = {2^j.Celh'^f/^/2 [21II2II23. It gives the approxi- 
mate size of the potential, that is, the wavefunction takes 
on its asymptotic scattering form for r 3> Xq. We have 
previously approximated the exact Born-Oppenheimer 
potential by the energy-dependent Fermi pseudopoten- 
tial Ullllil 

V,sir;E) = S{r)—r. (9) 

The dependence on collision energy E is due to the 
energy-dependent scattering length, defined as 



acff(-E) 



tan So{E) 



(10) 



where E = h^k^/2ii and So{E) is the s-wave collisional 
phase shift. S'-wave scattering predominates for ultracold 
collisions except for the case of identical fermions. 

The pseudopotential approximation is valid provided 
the van der Waals length scale Xq is less than the small- 
est harmonic oscillator width, xq ^ min{c?j^, dz} |2- We 
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assume that the energy shifts due to higher partial waves 
are neghgible in comparision to those of the s-wave. This 
tends to be true for ultracold colUsions because the cen- 
trifugal barrier heights are large compared to the collision 
energy. We also neglect inelastic losses, that is, the imag- 
inary part of the scattering length is much smaller than 
the real part 25]. 

Since spherical symmetry is broken by the anisotropic 
potential, does not commute with the Hamiltonian. 
Consequently, partial waves with the same parity with 
respect to I are coupled. While the projection of an- 
gular momentum on the z-axis does commute with the 
Hamiltonian, only its m; = eigenstates are affected by 
s-wave scattering. Thus we only compute the energies of 
even partial wave, rai — Q states. We do not consider odd 
partial waves, should they be present for distinguishable 
bosons, because they have negligible energy shifts in the 
limit of very low collision energy. 



III. METHOD OF EIGENVALUE SOLUTION 
AND SELF-CONSISTENT ENERGIES 

We need to solve the eigenvalue problem for the Hamil- 
tonian Eq. ^ self-consistently, because of the energy- 
dependent scattering length in the pseudopotential term. 
As in Ref. this is done in two steps. We first obtain 
the eigenvalues Ei{a/d, A)/fkj in scaled trap energy units 
for fixed values of A and the energy- independent scaled 
scattering length a/d. The self-consistent energy eigen- 
values for an actual system with acs{E) from Eq. H1U|) 
are then found graphically, for a given A, by superposing 
a plot of Ei{a/d,A) as a function of a/d and a plot of 
acs{E)/d, with E as the ordinate and Qos/d as the ab- 
scissa (cf. Fig. 1 and Fig 2 of 0)- The points where the 
curves intersect determine the self-consistent energies. 

In the rest of this section we focus on obtaining 
Ei{a/d,A), since aos{E) can be found from a standard 
free-space scattering calculation. We use the partial wave 
expansion of the wavefunction,and expand each partial 
wave, with the exception of s-waves, in isotropic har- 
monic oscillator eigenfunctions of frequency lo. For the 
s-wave part, we use the analytic eigenfunctions of the 
isotropic harmonic oscillator with a Fermi pseudopoten- 
tial proportional to scattering length a These au- 
tomatically incorporate the singular nature of the wave- 
function at the origin. Since m; = we can set the 
spherical coordinate (f) = 0, and consequently we write 

tp{r,0;a) = CnoQnjr; a)Yoo {9, 0) 

n 

+ ^^c„zi?„i(r)rjo(e,0) (11) 

n l>0 



Rniir) = 



r(n-f-;-|- f) 



(13) 

where the summations are over all nonnegative integers 
n and even I > 0. Here U is the (second) Kummer con- 

fluent hypergeometric function, L„ ^ are the Laguerre 
polynomials, and F is the Gamma function 26|. The 
nonintegral s-wave quantum numbers Vn for an isotropic 
trap are determined by 



a 
d 



1 , r(i/„ 

— tan TTVn — -. — 

2 V{Un 



1) 



2> 



(14) 



This equation also is used in the calculation of the deriva- 
tive in Eq. 

We need the Hamiltonian matrix elements in the basis 
used in Eq. lll|l . The isotropic part of the Hamiltonian 
Eq. Q gives only diagonal matrix elements: 



H, 



(0) 



(15) 



for / > 0. Only the s-wave diagonal matrix elements. 



(16) 



are affected at low energy by atom-atom interactions pro- 
portional to the scattering length. The anisotropic part 
of the Hamiltonian, Eq. I^J contributes both diagonal 
and off-diagonal matrix elements: 



H, 



= 

d\n'l' 



5 



(17) 



for all principal quantum numbers n, n' and even par- 
tial waves l^V . The derivation of these matrix elements 
is given in the Appendix. The diagonalization of the 
Hamiltonian matrix is straightforward with sparse matrix 
eigenvalue routines. For the most extreme anisotropics 
considered in this paper (as small as 0.01 and as large as 
100), we required a maximum / = 600 and a maximum 
n — 600 to compute the lowest few eigenvalues for all 
values of scattering length. We have checked that the 
correct solutions are approached as A — > 1. 



IV. QUASI-ONE-DIMENSIONAL TRAP 

We have found solutions in the quasi-one-dimensional 
regime ^ <C 1. Figure^shows Ei{a/d,A — 0.01) versus 
a/d. The figure also shows the eigenvalues El^{a/d) for 
a purely one-dimensional model, corresponding to inter- 
action via a delta function in z and trapping along z only. 
Following Ref. 16], the eigenvalues are 



Qn{r]a) 



2a 
7^ 



di^n , , f 3 



(12) 



3\5, A\ 
2 + 



l + A\ 2v. 



.ID 



(18) 



4 



FIG. 1: Energy eigenvalues versus scattering length at 
A = 0.01 from Fermi-pseudopotential (solid line) and one- 
dimensional scattering theory (dashed line).(j/ = tj/27r). 



FIG. 3: Energy eigenvalues versus scattering length at 
A = 100 from Fermi-pseudopotential (solid line) and two- 
dimensional scattering theory (dashed line).(j/ = tj/27r) 



FIG. 2; Wavefunction rxp{x,0,z) in Cartesian coordinates 
X — rsm6,y = 0, 2 = rcosO, corresponding to second lowest 
energy eigenvalue for A — 0.01 and a/d = —25. (b) is a close 
up of (a). All lengths are expressed in trap units of d and '0 
in units of d~^^^ . 



where the v}^ satisfy 



tan TTv, 



ID 



r(t./^ + 1) 



-.ID 



(19) 



The one-dimensional interaction, hujg^^6{z/d), is related 
to the three-dimensional scattering length through 



1 - 1.4603 (i)'(l + #)-- 



(20) 



Our eigenvalues agree well with those from the one- 
dimensional model for E > Q.Shui. To get this agreement, 
it is crucial to include the rcnormalization in the denom- 
inator of Eq. (pn|l . The levels for E > 1.23fta; represent 
trap levels aligned along the weak trapping direction with 
spacing about 2hujz = 0.024/iw. In Fig. there is no dif- 
ference within our numerical accuracy between our low- 
est three trap levels and those from the one-dimensional 
model maximum energy difference, while the difference 
for the highest trap level shown is 0.002huj. The eigen- 
values for large positive or negative scatttering length ap- 
proach the same asymptotic values. The lowest energy 
state is not predicted accurately by the one-dimensional 
model; we discuss this special state further in Sec. IVII 

The wavefunction ripi^x = rsin^, y = 0, z = rcosO) , 
corresponding to the second lowest energy, with A — 0.01 
and a/d= —25, is plotted in Fig. [5] In Fig.l^Ja) the one 
dimensional nature of the wavefunction is apparent on 
a scale large compared with d. The variation in x is 
approximately Gaussian, while scattering results in the 
dip along the line z = 0. A close-up of the wavefunc- 
tion in Fig. [2Ib) reveals how the scattering crosses over 
to a isotropic three-dimensional character at short inter- 
atomic distance. The function rijj is finite at the origin 
as a consequence of the pseudopotential scattering. 

We have compared some of our results to a recent dif- 
fusion quantum Monte-Carlo study on the ground state 
of interacting bosons in elongated traps ^2} ■ Our lowest 
trap state agrees with the ground state of that method to 
within one percent for anisotropics in the range 1 > A > 
0.01 for a fixed positive value of the scattering length, 
even though that study assumed a hard-core potential of 
size a. 



V. QUASI-TWO-DIMENSIONAL TRAP 

We have also computed solutions in the quasi-two- 
dimensional regime, 1. Figure O shows Ei{a/d,A — 
100) versus a/d. The figure also shows the eigenval- 
ues Ef^{a/d) for a purely two-dimensional model, cor- 
responding to interaction via zero-range two-dimensional 
scattering and trapping in the x, y-plane only. Again fol- 
lowing 16], the eigenvalues of a two-dimensional trap are 



A 

- + 2K 



2D 



huj . 



where the 



,2D 



satisfy 



^ » ) g2D 



(21) 



(22) 



and 5* is the digamma function 26] . The two- 
dimensional scattering is mediated by an interaction 
strength related to the three-dimensional scattering 
length through 



' 3 

\2-n'-^ 



A 



^(1 + ^ 



1 



' 3 

V 271^ 



\n{^)A\ 



1 



2 , 



(23) 



This expression for is derived by simple algebra from 
equations in Ref. |l9| . 

Our eigenvalues agree well with this model for E > 
O.Sfto;. To obtain this agreement it is crucial to include 
the rcnormalization in the denominator of Eq. I|23|l . 
The trap levels above E « 0.9fta; are spaced by about 
2hw±^ = 0.03462ftwj^. In Fig. El the difference between 
our lowest trap level and that from the two-dimensional 
model is 0.0015huj at a/d = ±20; the difference for the 
highest trap level shown is 0.002^01 at the same values of 
scattering length. The eigenvalues for large positive or 
negative scattering length approach the same asymptotic 
values. Again we note the exception that the lowest en- 
ergy state is not predicted well by the renormalized two- 
dimensional model (see Sec. IV1() . A feature of the two- 
dimensional physics seen in Fig.^lis that for \a\ <C d, all 
eigenvalue curves except the lowest have nearly the same 
slope. 

A sample wavefunction rip^x = r sin 0,7/ = 0,z = 
r cos 9) , corresponding to the second lowest energy, at 
A = 100 and a/d = -25, is plotted in Fig.|3| (Note that 
the X and z axes are interchanged from Fig. |21) The 
two-dimensional nature of the physics is immediately ap- 
parent on a scale large compared with d. The variation in 
z is approximately Gaussian, with a dip along x = 0. The 
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FIG. 4: Wavefunction rijj in Cartesian coordinates x — 
rsin6,y = 0,z = rcosO, corresponding to second lowest en- 
ergy eigenvalue for A — 100 and a/d = —25. (b) is a close up 
of (a). All units are as in Fig. |21 



FIG. 5: Wavefunction r'i/) in Cartesian coordinates x = 
rsin0,y — 0,z — rcosO, corresponding to lowest energy 
eigenvalue for A — 0.01 and a/d = —25. All units are as 
in Fig.H 

close-up Fig.^b) reveals how the scattering crosses over 
to a isotropic three-dimensional character at short inter- 
atomic distance, as in the quasi-onc-dimensional case. 

VI. TRAP-INDUCED BOUND STATE 

Curiously, for large magnitudes of the scattering 
length, the lowest energy eigenvalue appears near fiuj/2 
for all values of trap anisotropy. This behavior can be ex- 
plained by examining the wavefunctions of these states 
plotted in Figs. El and El at a = —25, for A = 0.01 and 
100 respectively. Both wavefunctions are more nearly 
isotropic than those corresponding to higher levels at the 
same parameters (see Figs. |2 and 0J, and their extent is 
roughly the mean trap length scale d. Similar wavefunc- 
tions are obtained for large positive scattering length. 

The energy of the lowest state is plotted in Fig. [71 for 
a/d — ±100 as a function of anisotropy. The shape of 
this curve can be understood from a perturbative picture 
in the anisotropic interaction H'-^\ defined in Eq. O. To 
zeroth-order in this picture the state is the lowest level of 
the isotropic trap of frequency uj and a pseudopotential 
with scattering length a. (Recall that for an isotropic 
trap with a — > ~oo the lowest energy is huj/2.) There 
is no first-order correction in H^^\ while to second-order 
the energy gets a small correction proportional to A^. 
Indeed, we find that a quadratic fit in A at a/d = —100 
results in Eo/huj = 0.5054 — 0.0562A^; compare with the 
zeroth-order energy 0.5056fiw for the lowest state in an 
isotropic trap at the same scattering length. Similarly, at 
a/d = 100 we find Eo/hw = 0.4941 - 0.0549A2 compared 
with the zeroth-order energy 0.4943/ia; in the isotropic 
case. 

When neglecting the pseudopotential, th e energy of t he 
lowest state is E{a = 0)/huj= {l+A/2)/y^2/3 + A'^/3 > 
\fi/2, so we should properly consider a state to be bound 
when its energy is lower than this. Thus we denote such a 
state the trap-induced bound state. One can also think of 
it as an artificial molecule with an extent given roughly by 



FIG. 6: Wavefunction r-f/; in Cartesian coordinates x — 
rsm9,y — 0,z — rcosO, corresponding to lowest energy 
eigenvalue for A = 100 and a/d = —25. All units are as 
in Fig.m 



FIG. 7: Lowest-energy eigenvalue versus trap anisotropy A 
(logarithmic scale) at a/d = —100 (triangles) and a/d — 100 
(squares). Fits up to first order in A^ for a/d — —100 (solid 
line) and a/d = 100 (dashed line) are shown; see text for 
coefficients of fits. The fits verify that the wavefunction is 
nearly isotropic for all A. [v = u/2ti) 



FIG. 8: Energy eigenvalues versus magnetic field B for two 
Na atoms in the lowest hyperfine level in a axisymmetric trap 
with V = u/2ti — 500kHz and anisotropy A = 0.01. The 
dashed line shows the energy of the lowest trap level when 
the interatomic interaction is neglected. 

the size of the tightest trap dimension. As the scattering 
length approaches zero from the positive side, we recover 
the usual molecular bound state with an energy below 
zero. 



VII. FESHBACH RESONANCE IN 
QUASI-ONE-DIMENSIONAL TRAP 

One possible way of varying the atomic interaction 
strengths experimentally is through the use of a tun- 
able Feshbach resonance state. Consequently we describe 
a quasi-one-dimensional magnetically tunable Feshbach 
resonance by using the self-consistent energy method 
with the eigenvalues of Fig.Q] We consider two Na atoms 
in their lowest hyperfine levels, for which an s-wave Fes- 
hbach resonance occurs near 90.9 mT [28L l29l l30l|. The 
scattering length is highly energy- and magnetic-field- 
dependent near the resonance. We use the effective scat- 
tering length from a close-coupling calculation, as de- 
scribed in Ref. 0, 1131 • The trap frequencies are taken 
to be 10^/211 = 612 kHz, lOz/2tt = 6.12 kHz, so the trap 
anisotropy is ^ = 0.01 and lo/2tt — 500 kHz. Using 
the procedure outlined at the beginning of Sec. IIIII we 
predict the eigenvalues as a function of applied magnetic 
field near the resonance in Fig. jSl As the magnetic field 
is tuned through resonance, the lowest state goes contin- 
uously from a molecular state with i? < 0, to the trap- 
induced bound state E hiu /2, to the lowest quasi-one- 
dimensional trap state at i? « y/3/2huj w l.23huj. For 
the trap frequencies chosen here one can change from a 
molecular bound state to a trap state by tuning the mag- 
netic field 0.01 mT. The trap states are smoothly shifted 
up by 2hujz as the magnetic field is increased. 



VIII. CONCLUSION 

We have argued that an energy-dependent pseudopo- 
tential approach may be used to calculate eigenvalues 
of two ultracold atoms colliding in an axisymmetric har- 
monic trap. Furthermore, we have numerically solved for 
the eigenvalues of the axisymmetric trap with an s-wave 
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pseudopotential interaction proportional to a scattering 
length. These results can be considered a generalization 
of the isotropic trap case previously solved 0, 13 • Our 
results show that one- and two- dimensional interaction 
regimes can be reached, but that the interactions become 
renormalized by the tight trapping potential when the 
magnitude of the effective scattering length is large com- 
pared with a mean trap length. Remarkably, in the case 
of scattering length of large magnitude, we find a nearly 
isotropic state with energy near hio /2 for all values of trap 
anisotropy. This is a trap-induced bound state. The size 
of state is controllable by the tightest trap frequency. We 
show by an example that this state can be reached with 
the current techniques of magnetically tunable Feshbach 
resonance applied to atoms in an optical lattice. 

The numerical techniques used in this article may also 
be useful when the interatomic interaction becomes suffi- 
ciently long-range compared with the tightest trap direc- 
tion or is anisotropic such that the s-wave pseudopoten- 
tial approach becomes insufficient. This is of im por tance, 
for example, for dipole-dipole interactions 0jl2H23l- The 
possibility of manipulating shape resonances (such as the 
d-wave resonance in Na collisions) with the trap should 
not be overlooked. The effective-scattering-length model 
and partial wave expansion in our numerical technique 
should also solve the eigenvalue problem for atoms collid- 
ing in separated traps, as can occur in a state-dependent 
optical lattice. This is particularly important for the pro- 
posals on quantum computing with neutral atoms. We 
are continuing work in this area. 
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(2Z;00|2Z;Z'0) shows that 



hi 



(2Z-3)(2Z + 3) 



1+2,1 



3(/ + l)(Z + 2) 



2(2/ -I- 3)V(2/ + l)(2/-f 5) 



(25) 



(26) 



while all other angular matrix elements are zero. 

This leaves the radial factor of the matrix element to 
be computed. Throughout the remainder of the appendix 
kets refer to the radial part of the basis functions only, 
so that for (cf. Eqs. (O and lO) 



{r\nl) 
(?inO) 



Rni{r) 
Q„{r;a) . 



In this notation 



H 



(1) 

nl\7i'l' 



^jiU}^Mw{niy\n'l') 



(27) 
(28) 



(29) 



The analytic evaluation of the radial matrix elements 
{nl\r'^\n'l) is most conveniently carried out using n- and l- 
ladder operators. From the factorization method applied 
to the radial Schrodinger equation 33] , the n-raising and 
n-lowering operators are 



(30) 



b-i\nl) = yj2n{2n + 2l + \)\n-\,l) (31) 
= \/2n{2n + 2l + l)\nl), (32) 



while the Z-raising operator is 



such that 



or r 



Cl\nl) = -2V^|n-l,; + l). 



(33) 



(34) 



Appendix: Matrix elements of anisotropic potential 

In this appendix we evaluate the matrix elements of the 
anisotropic potential term in the basis of partial 

waves, isotropic harmonic oscillator functions and the ir- 
regular s-wave oscillator eigenfunctions. 

The partial wave expansion applied to the angle- 
dependent factor of the Hamiltonian is evaluated with 
the three-spherical-harmonic formula. 



21 + 1 



(2?;00|2/;/'0) 



(24) 



\ 21' + 1 

Evaluation of the Clebsch-Gordan coefficient 



(Note that these operators are applied to the normalized 
radial eigenfunctions.) 

For I = I' > Q, we use the fact that 



3 1 - 1 - 
= 2n -f Z H 6+, , , 6", 



and orthonormality to find 

{nl\r^\nl) = 2n + I + ^ 
{n + l,l\r'^\nl) = {nl\r^ \n + 1 , 1) 



(35) 



(36) 



= --V2(n + l)(2n + 2/ + 3), (37) 

with all other equal-Z matrix elements vanishing. (We 
recognize Eq. H36|) as a consequence of the quantum virial 
theorem for the isotropic harmonic oscillator.) 



7 



For I' = 1 + 2 (but ? 7^ 0), we begin with 

,, ,,,.„ (,.,|.^£^.£^|„- + 2.0 
4V('i' + l)K + 2) 

and use the operator identity 

r2£+,£+ = \-b-^,, + 2{n + l) + 2 



■'n+2 



j + 2{n + 2) 



to obtain the non-vanishing matrix elements 



1 



(38) 



(39) 



{nl\r^\n,l + 2) = -v/(2n + 21 + 3)(2n + 2Z + 5)(40) 



(n + l,Z|r2|n,/ + 2) = - v/2(n + l)(2n + 2/ + 5X41) 
(n + 2,Z|r2|n,Z + 2) = + + 2). (42) 



For the special case of Z = 0, we use the expansion of 
the irregular solutions in terms of \nl) and apply the 
above matrix elements. This results in 



/2r(n' + |) /5r/„ 



{n0\r^\n'2) c, 
y ni [n' + 1) \ oa 

1 2 1 



n' — Vn n' + \ — Vn n' + 2 — v^, 



(43) 



(The matrix elements (nO|r^|n'0) are not needed since 
/oo = 0.) 
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